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We investigate the conditions under which cosmological variations in physical 'constants' and 
scalar fields are detectable on the surface of local gravitationally-bound systems, such as planets, in 
non-spherically symmetric background spacetimes. The method of matched asymptotic expansions 
is used to deal with the large range of length scales that appear in the problem. We derive a sufficient 
condition for the local time variation of the scalar fields driving variations in 'constants' to track their 
large-scale cosmological variation and show that this is consistent with our earlier conjecture derived 
from the spherically symmetric problem. We perform our analysis with spacetime backgrounds 
that are of Szekeres-Szafron type. They are approximately Schwarzschild in some locality and 
free of gravitational waves everywhere. At large distances, we assume that the spacetime matches 
smoothly onto a Friedmann background universe. We conclude that, independent of the details of 
the scalar-field theory describing the varying 'constant', the condition for its cosmological variations 
to be measured locally is almost always satisfied in physically realistic situations. The very small 
differences expected to be observed between different scales are quantified. This strengthens the 
proof given in our previous paper that local experiments see global variations by dropping the 
requirement of exact spherical symmetry. It provides a rigorous justification for using terrestrial 
experiments and solar system observations to constraint or detect any cosmological time variations 
in the traditional 'constants' of Nature in the case where non-spherical inhomogeneities exist. 

PACS Nos: 98.80.Es, 98.80.Bp, 98.80.Cq 

I. INTRODUCTION 



Over the past few years there has been a resurgence of observational and theoretical interest in the possibility that 
some of the fundamental 'constants' of Nature might be varying over cosmological timescales [1]. In respect of two 
such 'constants', the fine structure constant, a, and Newton's 'constant' of gravitation, G, the idea of such variations 
is not new, and was proposed by authors such as Milne [2] , Dirac [3] , and Gamow [4] as a solution to some perceived 
cosmological problems of the day [5]. At first, theoretical attempts to model such variations in constants were rather 
crude and equations derived under the assumption that constants like G and a are true constants were simply altered 
by writing-in an explicit time variation. This approach was first superseded in the case of varying G by the creation 
of scalar-tensor theories of gravity [6], culminating in the standard form of Brans and Dicke [7] in which G varies 
through a dynamical scalar field which conserves energy and momentum and contributes to the curvature of spacetime 
by a means of a set of generalised gravitational field equations. More recently, such self-consistent descriptions of 
the spacetime variation of other constants, like a [8, 9], the electroweak couplings [10], and the electron-proton mass 
ratio, /x, [11] have been formulated although most observational constraints in the literature are imposed by simply 
making constants into variables in formulae derived under the assumption that are constant. 

The resurgence of interest in possible time variations in a and [i has been brought about by significant progress 
in high-precision quasar spectroscopy. In addition to quasar spectra, we also have available a growing number of 
laboratory, geochemical, and astronomical observations with which to constrain any local changes in the values 
of these constants [12]. Studies of the variation of other constants, such as G, the electron-proton mass ratio, 
(i = m e /m pri and other standard model couplings, are confronted with an array of other data sources. The central 
question which this series of papers addresses is how to these disparate observations, made over vastly differing scales, 
can be combined to give reliable constraints on the allowed global variations of a and the other constants. If a varies 
on cosmological scales that are gravitationally unbound and participate in the Hubble expansion of the universe, 
will we see any trace of this variation in a laboratory experiment on Earth? After all, we would not expect to 
find the expansion of the universe revealed by any local expansion of the Earth. In Paper I [13], we examined this 
question in detail for spherically symmetric inhomogeneous universes that model the situation of a planet or a galaxy 
in an expanding Friedmann-Robertson- Walker (FRW)-like universe. In this paper we relax the strong assumption 
of spherical symmetry and examine the situation of local observations in a universe that contains non-spherically 
symmetric inhomogeneity. Specifically, we use the inhomogeneous metrics found by Szekeres to describe a non- 
spherically symmetric universe containing a static star or planet. As in Paper I, we are interested in determining the 
difference (if any) between variations of a supposed 'constant' or associated scalar field when observed locally, on the 
surface of the planet or star, and on cosmological scales. 
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When a 'constant', C, is made dynamical we can allow it to vary by making it a function of a new scalar field, 
C — > C(4>), that depends on spacetime coordinates: cp = <t>(,%>t). It has become general practice to combine take all 
observational bounds on the allowed variations of C((f>). This practice assumes implicitly that any time variation of C, 
on or near the Earth, is comparable to any cosmological variation that it might experience, that is to high precision 

<j>{x,t) ~j>c(t), (1) 

for almost all locations x, where 4> c is the cosmological value of <j>. This assumption is always made without proof, 
and there is certainly no a priori reason why it should be valid. Strictly, 4> mediates a new or 'fifth' force of Nature. 
If the assumed behaviour is correct then this force is unique amongst the fundamental forces in that its value locally 
reflects its cosmological variation directly. 

In this series of papers we are primarily interested in theories where the scalar field, or 'dilaton' as we shall refer to 
it, <f), evolves according to the conservation equation 

□0 = B (0)KT-y (0), 

where T is the trace of the energy momentum tensor, T = (with the contribution from any cosmological constant 
neglected). We absorb any dilaton-to-cosmological constant coupling into the definition of V(<j>). The dilaton-to- 
matter coupling B((f>) and the self-interaction potential, V(<p), are arbitrary functions of 4> and units are defined by 
k = 8ttG and c = h = 1. This covers a wide range of theories which describe the spacetime variation of 'constants' 
of Nature; it includes Einstein-frame Brans-Dicke (BD) and all other, single-field, scalar-tensor theories of gravity 
[7, 9, 14, 15]. In cosmologies that are composed of perfect fluids and a cosmological constant, it will also contain 
the Bekenstein-Sandvik-Barrow-Magueijo (BSBM) theory of varying a, [9], and other single-dilaton theories which 
describe the variation of standard model couplings, [15]. We considered some other possible generalisations in [13]. It 
should be noted that our analysis and results apply equally well to any theory which involves weakly-coupled, 'light', 
scalar fields, and not just those that describe variations of the standard constants of physics. 

In first paper of this series, [13], we determined the conditions under which condition 1 would hold near the surface 
of a virialised over-density of matter, such as a galaxy or star, or a planet, such as the Earth, under the assumption 
of spherical symmetry. We chose to refer to this object as our 'star'. In Paper I, matched asymptotic expansions 
were employed to analyse the most general, spherically- symmetric, dust plus cosmological constant embeddings of the 
'star' into an expanding, asymptotically homogeneous and isotropic spherically symmetric universe. We proved that, 
independent of the details of the scalar-field theory describing the varying 'constant', that 1 is almost always satisfied 
under physically realistic conditions. The latter condition was quantified in terms of an integral over sources that can 
be evaluated explicitly for any local spherical object. 

In this paper we extend that analysis, and our main result, to a class of embeddings into cosmological background 
universes that possess no Killing vectors i.e. no symmetries. The mathematical machinery that we use to do this is, 
as before, the method of matched asymptotic expansions, employed in [13], where the technical machinery is described 
in detail. A summary of the results obtained there can also be found in [16]. 

This paper is organised as follows: We shall firstly provide a very brief summary of the method of matched 
asymptotic expansions used here. In section II we will introduce the geometrical set-up that we will use. We will 
be working in spacetime backgrounds of Szekeres-Szafron type [17, 18]. We describe theses particular solutions of 
Einstein's equations briefly in section II and then in greater detail in section III. In section IV we extend the analysis 
of [13] to include non-spherically symmetric backgrounds of Szekeres-Szafron type. In section V, we consider the 
validity of the approximations used, and state the conditions under which they should be expected to hold. In section 
VI we perform the matching procedure (as outlined below), and extend the main result of [13] to Szekeres-Szafron 
spacetimes. We consider possible generalisations of our result in section VII before considering the implications of 
the results in the section VIII. 

We will employ the method of matched asymptotic expansions [19, 20]. We solve the dilaton conservation equations 
as an asymptotic series in a small parameter, 6, about a FRW background and the Schwarzschild metric which 
surrounds our star. The deviations from these metrics are introduced perturbatively. The former solution is called 
the exterior expansion of 4>, and the latter the interior expansion of <j>. The exterior expansion is found by assuming 
that the length and time scales involved are of the order of some intrinsic exterior length scale, Le- Similarly in the 
interior expansion we assume all length and time scales we be of the of L/, the interior length scale. Neither of the two 
different expansions will be valid in both regions. In general, we define 8 := Lj/Le <C 1. This means that in general 
only a subset of our boundary conditions we will be enforceable for each expansion, and as a result both the interior 
and exterior solutions will feature unknown constants of integration. To remove this ambiguity, and fully determine 
both expansions, we used the formal matching procedure. The idea is to assume that both expansions are valid in 
some intermediate region, where length scales go like Li nt — LfL 1 ^ 01 , for some a e (0, 1). Then by the uniqueness 
property of asymptotic expansions, both solutions must be equal in that intermediate region. This allows us to set 
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the value of constants of integration, and effectively apply all the boundary conditions to both expansions. A fuller 
discussion of this method, with examples, and its application in general relativity is given in [13]. 

II. GEOMETRICAL SET-UP 

We shall consider a similar geometrical set-up to that of Paper I. We assume that the dilaton field is only weakly 
coupled to gravity, and so its energy density has a negligible effect on the expansion of the background universe. 
This allows us to consider the dilaton evolution on a fixed background spacetime. We will require this background 
spacetime to have the same properties as in Paper I, but with the requirement of spherical symmetry removed: 

• The metric is approximately Schwarzschild, with mass m, inside some closed region of spacetime outside a 
surface at r = R s . The metric for r < R s is left unspecified. 

• Asymptotically, the metric must approach FRW and the whole spacetime should tend to the FRW metric in 
the limit m —* 0. 

• When the local inhomogeneous energy density of asymptotically FRW spacetime tends to zero, the spacetime 
metric exterior to r = R s must tend to a Schwarzschild metric with mass m . 

We will also limit ourselves to considering spacetimes in which the background matter density satisfies a physically 
realistic equation of state, specifically that of pressureless dust (p = 0). We also allow for the inclusion of a cosmological 
constant, A. The set of all non-spherical spacetimes that satisfy these conditions is too large and complicated for us 
to examine fully here; and such an analysis is beyond the scope of this paper. We can simplify our analysis greatly, 
however, we specify four further requirements: 

1. The flow-lines of the background matter are geodesic and non-rotating. This implies that the flow- lines are 
orthogonal to a family of spacelike hypersurfaces, S t . 

2. Each of the surfaces St is conformally flat. 

3. The Ricci tensor for the hypersurfaces St, ^R a b, has two equal eigenvalues. 

4. The shear tensor, as defined for the pressureless dust background, has two equal eigenvalues. 

The last three of these conditions seem rather artificial; however, when the deviations from spherical symmetry are 
in some sense 'small' we might expect them to hold as a result of the first condition. In the spherically symmetric 
case, condition 1 implies conditions 2, 3 and 4. In the absence of spherical symmetry, these conditions require the 
background spacetime to be of Szekeres-Szafron type, containing pressureless matter and (possibly) a cosmological 
constant. The conditions (1-4) combined with the background matter being of perfect fluid type provide an invariant 
definition of the Szekeres-Szafron class of metrics that is due to Szafron and Collins [21, 22]. 

We have demanded that the 'local' or interior region be approximately Schwarzschild. The intrinsic length scale of 
the interior is defined by the curvature invariant there: 

R^ ^ 

L I = (l2 RabcdR^) = S 1/2 - (2) 

(2m) ' 

The exterior (or cosmological) region is approximately FRW, and so its intrinsic length scale is proportional to the 
inverse square root of the local energy density: 1/y/ne + A, where e is the matter density. In accord with current 
astronomical observations, we assume that this FRW region is approximately flat, and so we set our exterior length 
scale appropriate for the present epoch, t = t a , by the inverse Hubble parameter at that time: 

Le = 1/Ho- 

We can now define a small parameter by the ratio of the interior and exterior length scales: 



S = Li/L E - 
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III. SZEKERES-SZAFRON BACKGROUNDS 



In 1975 Szekeres [17] solved the Einstein equations with perfect fluid source by assuming a metric of the form: 

ds 2 = dt 2 - e 2a dr 2 - e 2 ' 3 {dx 2 + dy 2 ) , 

with a and (3 being functions of (t,r,x,y). The coordinates where assumed to be comoving so that the fluid-flow 
vector is of the form: m m = 5^ ; This implies p = p(t) and the acceleration ii^ = 0. Szekeres assumed a dust source 
with no cosmological constant, p = 0, although his results were later generalised to arbitrary p(t) by Szafron [18] 
and the explicit dust plus A solutions were found by Barrow and Stein-Schabes [24]. In general, these metrics have 
no Killing symmetries [23]. Spherically-symmetric solutions of this type with a(r, t) and f3(r,t) were, in fact, first 
discussed by Lemaitre [25] and arc usually referred to as the Tolman-Bondi models [26] ; much of the analysis of Paper 
I assumed a Tolman-Bondi background. 

The Szekeres-Szafron models can be divided into two classes: /3 r = and /3 r ^ 0. Both classes include all FRW 
models in their homogeneous and isotropic limit; however, only the latter 'quasi-spherical' class includes the external 
Schwarzschild solution. Since we want to have some part of our spacetime look Schwarzschild we will only consider 
the /? r quasi-spherical solutions. We will also limit ourselves to spacetimes with a cosmological constant, [24], 
so in effect the total pressure is p = —A. These universes contain no gravitational radiation as can be deduced from 
the existence of Schwarzschild as a special case which ensures a smooth matching to Schwarzschild, which contains 
no gravitational radiation. With these restrictions, a and (5 are given by: 

e = <t>(t,r)e Hr < x < v \ (3) 
e a = h(r) e -^ r ' x ' v) (e ) r , (4) 



where r) satisfies: 



e~ v = A(r)(x 2 +y 2 ) + 2B 1 (r)x + 2B 2 (r)y + C(r), (5) 



$ 2 = -k(r) + 2M(r)/$ + ^A$ 2 . 



The functions A{r), -Bi(r), B 2 (r), C(r), M(r), k(r) and h(r) are arbitrary up to the relations: 

\E(r) := AC B\ B\ = \ [^ 2 (r) + k(r) . 

The surfaces (t, r) — const have constant curvature E(r). We will require that the inhomogeneous region of our 
spacetime is localised, so that it is by some measure finite. This implies that the surfaces of constant curvature must 
be closed; we must therefore restrict ourselves to only considering backgrounds where E > 0. Whenever this is the 
case, we always can rescale the arbitrary functions so that E can be set equal to 1 by the rescalings 

A(r) := A{r)/^E(rj, B 1 (r) := B^/^W), B 2 (r) := B 2 (r)/^E{r), C(r) := C{r)/y/E(rj, e v := ^Ee" 
k := k(r)/E(r), h~ 2 := h{r)- 2 /E{r) = 1 - k(r), R(t, r) := $(r, t)/VE, M(r) = M(r)/E 3/2 . 

These transformations can be viewed as the 'gauge-fixing' of arbitrary functions. In this gauge, R(t, r) is a 'physical' 
radial coordinate, i.e. the surfaces (i, r) = const have surface area AttR 2 and the metric becomes 

„ „ (1 + v rR) 2 R 2 dr 2 „ „ . „ „, 

ds 2 = dt 2 - f ' - r R 2 e 2v (dx 2 + dy 2 ) , 

1 - k(r) v ' 

where e~ v = A(r)(x 2 + y 2 ) + 2B 1 (r)x + 2B 2 (r)y + C(r) and AC - B\ - B\ = ±, and v. R := u tr /R, r and: 

R 2 t = -k(r) + 2M(r)/R+ -AR 2 . 

In this quasi-sphcrically symmetric subcase of the Szekeres-Szafron spacetimes the surfaces of constant curvature, 
(t,r) = const, are 2-sphcres [27]; however, they are not necessarily concentric. In the limit v r — > 0, the (t, r) = const 
spheres becomes concentric (see fig. 1). We can make one further coordinate transformation so that the metric on 
the surfaces of constant curvature, (t, r) = const, is the canonical metric on S 2 i.e. d# 2 + sin 2 9d(f> 2 : 

x^X = 2(A(r)x + B 1 (r)), 
y^Y = 2(A(r)y + B 2 (r)), 
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where X + iY = cot 0/2. This yields 

-tvky = x 2 +Y 2 + 1 = A z (r)cos0 + A x (r)sin0cos<£ + A y (r) sin sin ip, 

where we have defined: 

A,(r):=^, A,(r):=(^)'A A,(r) := (^)'a (6) 
With this choice of coordinates, the local energy density of the dust separates uniquely into a spherical symmetric 




< ■ f — : ; surfaces of 

constant curvature 
/ / .' (t,r)=const. 



cosmological region 



FIG. 1: The surfaces (t, r) are spheres, which are concentric to leading order in 8 in both the exterior and interior limits, 
part, e s , and and a non-spherical part, e ns : 

e = e s (t,R) +e ns (t,R,9,<p), 



where: 

c 



K£ s - -5*-, (7) 



RKR (2M , 

K£ ™ = ~ TTT. E 1 "55- ■ ( 8 ) 



l + v, R R \R 3 



.r 



We define M t R = M, r /R yr . Following the conventions of our previous paper we write define 

M := m + Z(r), 

where m is the gravitational mass of our 'star'. 

A. Exterior Expansion 

As a result of the way that the inhomogencity is introduced in these models, we want the FRW limit to be 'natural' 
, that is for the 0(3) orbits to become concentric in this limit; we therefore require v r ~ o(l) as 5 — > in the exterior. 
This follows from the requirement that the whole spacetime should become homogeneous in a smooth fashion in 
the limit where the mass of our 'star' vanishes: m — > 0. Put another way, the introduction of our star is the only 
thing responsible for making the surfaces of constant curvature non-concentric. We define, as in the previous paper, 
dimcnsionless 'radial' and time coordinates appropriate for the exterior by 



r = H t, p = H r. 
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The exterior limit is defined by 5 — > with r and p fixed. In the exterior region we find asymptotic expansions in 
this limit. According to our prescription, we write 

H Z(p) ~ Ul mP 3 + 5P Zl (r) + o(6 p ), 

and 

H - 1 A i ~5'i i (p) + ?(5 a ). 

Since H 2 (^-) R ~ £), we have that: H^ 2 K e ns ~ 0((5p+ s , (5 1+;s ) whereas #~ 2 k£ s - £). Thus, 

the non-spherical perturbation to the energy density is always of subleading order compared to the first order in 
spherical perturbation. The first-order, non-spherical, metric perturbation appears at 0(S S ); however, since this is 
equivalent to a coordinate transform on (r, 0, <p) and the dilaton field, 0, is homogeneous to leading order in the 
exterior, this perturbation does not make any corrections to the dilaton conservation equation at 0(5 S ). Thus, both 
at leading order, and at next-to-leading order, both the energy density and the dilaton field will behave in the same 
way as in the spherically-symmetric Tolman-Bondi case - with the possible addition of a non-spherically symmetric 
vacuum perturbation to the dilaton, </>, i.e. <fi <~ <f> s + 4> ns + o(S p ) where cj) s is the spherically symmetric solution and 
OFRw4>ns = 0. As in our previous paper, however, we are not especially interested in the exterior solution for <fi 
beyond zeroth order, just the effect of any background variation in <f> on what is measured on the surface of a local 
'star'. 



B. Interior Expansion 

We define dimcnsionlcss coordinates for the interior in the same way as we did for the spherically symmetric case: 

T = L-J l (t - t ) and £ = R/R a . 

We define the interior limit to be 8 — > with T and £ fixed, and perform out interior asymptotic expansions in 
this limit. To lowest order in the interior region, we write Z <~ 8 q R s fi\, and Aj := S q >Rj 1 b i , where i = {x,y, z}. 
The condition that ne > everywhere requires q 1 > q and then, to next-to-leading order, the interior expansion of 
4> will be the same as it was in the spherically-symmetric Tolman-Bondi case. We can potentially include a non- 
spherical vacuum component for <fi; however, this will be entirely determined by a boundary condition on R = R s 
and the need that it should vanish for large R. To find the leading-order behaviour of the 4>^t we need to know 
(j) at next-to-leading order. The only new case we need to consider therefore is when q' = q, In 
the spherically symmetric case we considered two distinct subclasses of the Tolman-Bondi models: the flat, k = 0, 
Gautreau- Tolman-Bondi spacetimes, [22, 28] and the non-flat, k ^ 0, Tolman-Bondi models with a simultaneous 
initial singularity. In Gautreau- Tolman-Bondi models the initial singularity is non-simultaneous from the point of 
view of geodesic observers. The latter class is the more realistic, since in the former the world-lines of matter particles 
stream out of the surface of our star at R — R s \r=r s i.e. i? t > 0, whereas in the simultaneous big-bang models we 
can demand that matter particles fall onto this surface i.e. R,t\R=R e < 0. With this choice, and if R s = 2m, the 
non-flat models properly describe the embedding of a black hole into an expanding universe, whereas the Gautreau- 
Tolman-Bondi model technically describes the embedding of a white-hole in the same universe. In this paper we 
shall, therefore, only give the results explicitly for the non-flat case - however, we can present a simple procedure to 
transform our results to the flat Gautreau case. 

We define 

n = (e 12 - 3T/2) 2/3 ; R s n = r + 0{8\ S 2 / 3 ). 

From the exact solutions we find: 

k(n) - S 2 / 3 k (1 + a^ifa) + o (<5 9 )) + O (<5 5 / 3 ) , 

where 

, , rm , 2m f it \ 2/3 
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We can remove the 0{5 2 ^) metric perturbation by a redefinition of the T coordinate, T — > T*: 



2m f -I _ -I _ ( 5 2 / 3 tt£> \ 

, /■« V V V \H t +8T) J 

Vl-SVW*=T+ / V ' ^d£'. 

To leading order we see that T ~ T* . The interior expansion of the metric, for q' = <j, is written: 

d*L - Rl (if^O + ^ )8 «. X) + ^"(e, X) + dx a dx b + o(Si). 

where and are given by: 

j^dx a dx b = T^dT* 2 - (d£ + r 1/2 dT*) 2 - £ 2 {d# 2 + sin 2 9dp 2 }, (9) 
j^dzW = -^d£dT* - ^M dT *2. (10) 

These are the same as in the spherically symmetric case. The non-sphcrically symmetric perturbation is given by 

j { ab nS dx a dy b = 2(b z cos9 + b x cos ip sin 6 + b y sin p sin 9)£dr) 2 (11) 
— 2£ 2 [b z sin 9 + b x cos <^(1 — cos 9) + b y sin p>(l — cos 0)] d6*d?7 
— 2£ 2 (1 — cos#) sin9(b x simp — b y cosp)dpdrj. 

The spherically symmetric part of the local energy density ne s is the same as it was in the Tolman-Bondi cases: 

s s " R s e /2 V 1/2 ' 

The non-spherically symmetric part is: 

2 6m (6 2 cos 9 + b x sin cos <f> + b y sin sin </>) 

Tl s K£ ns = - 



R s ^3/2^1/2 



and to ensure that the energy density is everywhere positive we need (jffl > 36^. 

IV. EXTENSION TO QUASI-SPHERICAL SITUATIONS 
A. Boundary Conditions 

We demand the same boundary conditions as before: as the physical radius tends to infinity, R — > oo, we demand 
that the dilaton tends to its homogeneous cosmological value: (f>(R,t) — > <fi c (t). This can be applied to the exterior 
approximation. In the interior, we demand that the dilaton-flux passing out from the surface of our 'star' at R = R s 
is, at leading order, parametrised by: 

-Rs (l - |r) d^ \^ Rs = 2mF (fa) = J ' dR'R ,2 B^(Mi'))^(R'), (12) 

where (f>o — <fra(R = R s )- The function F(<f>) can be found by solving the dilaton field equations to leading order 
in the R < R s region. If the interior region is a black-hole (R s = 2m) then we must have F{4>) — 0; otherwise we 
expect F{4>) ~ B^(<p). Without considering the sub-leading order dilaton evolution inside our 'star', i.e. at R < i? s , 
we cannot rigorously specify any boundary conditions beyond leading order. Despite this, we can guess at a general 
boundary condition by perturbing eq. (12): 

-i? 2 (l - ^] d R S{<P) _ = - 5 (V=99 RR ) d R <f>o _ + 2S(M)F (fa) + 2mF^(fa)S (fa) + smaller terms, (13) 
\ j^s / £ — Rs R — Rs 

where 6(X) is the first sub-leading order term in the interior expansion of X; M is the total mass contained inside 
£ < R s and is found by requiring the conservation of energy; and at t = to we have M — m. Only 5 (fa) remains 
unknown; however, we shall assume it to be the same order as 5((j>) and see that this unknown term is usually 
suppressed by a factor of 2m/R s relative to the other terms in eq. (13). 
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B. Interior Expansion 



In the spherically symmetric case we found that 



+8 q <t>} ' +o(5 q ). In the non-spherical case, where q' = q, 



we relabel 



a(i) 



and we have additional non-spherical modes: 



4 0) (£. T ) + S q (f> { i )s {L T) + 8 q 4> [ p z {£„ T) cos 9 + 5 q 4>f )x % T) sin 9 cos <p + 8 q 4>f )v (£, T) sin 9 sin if + o(S q ) 



where: 



+ 



+ 




(14) 



We can solve this order by order in 2m/R s , and to lowest order we find: 



2m 



b,* (0?) e 



\ Mi) 



2m 

R.s 



B,* (0?) 



WW) 



(15) 



9m 1 f <7- 

+ 1T F $>) 71 / d'TWtf' + + A£ + 0((2m/i2 t 



) 2 ) 



Since we arc interested in finding when and where the local time variation of <f> deviates from its cosmological value, 
we are chiefly concerned with the case q < 1. The matching condition then requires that we fix Di so that in the 
intermediate limit we have <j^p l <~ £™ with n < 1. The value of C, should be set by a boundary condition on R = R s . 
We cannot specify Cj exactly without further information about the interior of our 'star' in R < R s . If we assume 
that the prescription for the sub-leading order boundary condition given above is correct then we find: 



2m hi 



Ci = 



R s v 1/2 

m 
R. 



0((2m/R s ) 2 ) 



£=i 



■B^ ($) 



From now onwards we set Cj = 0, for simplicity; even when this is not correct we do not expect the magnitude of Ci 
or Ci_T to be larger than any of the other terms in <j^p l or (jffiji , respectively. The time-derivative of <j)"p l for fixed R 
is: 



} I,T 



4m 
Rs 



B, 



Q « £'5/2 



2m 



-5. 



£1/2 



1 

e 



£1/2 



D,tZ (16) 



In the next section we shall discuss what we require of the &, for the matching procedure to be valid. In section VI 
we will then use the matching conditions to find D and D t- 

We could also relax the requirement that the leading-order mode in 4> be spherically symmetric. At next-to-leading 
order these new modes would generate extra terms in (ffi . In general, an /-pole at leading order becomes an I + 1-polc 
at next-to-leading order. The magnitude of the extra time-dependence that is picked up is, however, the same each 
time. Hence, we restrict ourselves by taking the leading-order mode to be spherically symmetric for the time being. 
Note also that we can pass from the simultaneous big-bang case, to the spatially flat, 'Gautreau', case by setting 

k = and making the transform r/ — > \ = (£ 3 ^ 2 + 3T/2) 2//3 . This will also mean that 4>j^t - 



01, T- 



V. VALIDITY OF APPROXIMATIONS 



All of the conditions found in Paper I for the matching of the spherically symmetric parts of <f> to be possible 
still apply here. However, we must now satisfy some extra conditions that come from the requirement that the 
non-spherical parts should also be matchable. 

We assume that 6j (77) oc r\ di as r\ — > 00 for some di > 0. At order 8 q , the growing mode in the non-spherically 
symmetric part of the interior approximation will then grow like 8 q rj di+1 /£. In the intermediate, or matching, region 
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we have that 77, £ ~ 6 "for some a £ (0, 1). We require (j>i to have a valid asymptotic expansion this region. This 
implies that there exists some a £ (0, 1) such that, for each i, we have a — q/di > 0. 

In the exterior we shall write H^ 1 ^ ~ d Pi k(p) 1 where p[ > comes from the requirement that the 2-spheres of 
constant curvature become concentric in the exterior limit. As p — > we assume that h(p) oc p~^. We previously 
stated that Z <~ ^fl m p 3 + <5 p ^i + o(5 p ) in the exterior. We assume that as p — > 0, we have zi oc p~ m . Although we 
did not explicitly consider the exterior expansion of <j> we can now examine the behaviour of the leading-order non- 
spherically symmetric mode in the intermediate limit of that exterior expansion. We noted above that there will be 

no 0(5 Pi ) correction resulting from the Zj. The leading-order mode will therefore either go like max^ I 8 P+Vi Z\(p)li{p) 



if p < 1 or maxi (s 1+Pi (p)li(p)j otherwise, and p ~ ©((5 1 Q ) in the intermediate region. For the exterior expansion 
to be valid in the intermediate region we therefore require 

max(p- + (1 - a)(fi + to)) > —p if p < 1, 

i 

max(p- + (1 - a)fi) > -1 if p > 1. 

i 

These conditions on a are equivalent to the following: there exists a such that the interior expansion of R 2 ne ns is 
o(l) as S — > for all < a' < a where £, T ~ 0(<5 _a ), and the exterior expansion of R 2 ne ns is also o(l) as <5 — > for 
all < a" < a where p,r — r ~ 0((5 1_0! ). This suggests that the condition for the matching procedure to work, as 
far as the spherically non-symmetric modes are concerned, is simply that 

R 2 K£ ns -C 1 everywhere. 

We can also rephrase and generalise the conditions for the matching procedure to be possible w.r.t. the spherically 
symmetric modes (as found in [13]) in a similar fashion: for all a e (0, 1), and keeping LJL x ^ a {t — to), LJL l F T a R 
fixed, we have lim^o i? 2 «:Ae s = o(l) and lim^o 2(to + Z)/R = o(l). We can combine our two conditions by simply 
replacing Ae^ by Ae in the above expression. Strictly speaking, since a £ (0,1) (as opposed to [0,1), (0,1] or 
[0, 1]) we can also replace Ae by just e since R 2 k£frw is small everywhere outside the exterior region. For Szekeres 
backgrounds the first of these conditions implies the second everywhere outside the interior region. Therefore, the 
matching procedure is certainly possible to zeroth order, if: 

Va e (0, 1) : lim (R 2 ne(R,t)) = o(l) and lim (M(R,t)/R) = o(l) with {L a I L 1 E a {t - t ), LfL 1 E a R} fixed, 

where M(R,t) is the gravitational mass inside the surface (t, R) — const. Equivalently, in any intermediate region 
the background spacetime is asymptotically Minkowski as 5 — > 0: everywhere which is not in either the interior or 
exterior regions can be considered to be a weak-field perturbation of Minkowski spacetime. The power of our method 
is that we do not require this to be true of the interior and exterior regions. So long as this condition holds in the 
intermediate region, we can match the zeroth-order approximations in some region and find the circumstances under 
which condition (1) holds by comparing the relative sizes of the derivatives <j> c . t and 4>f\. 

VI. MATCHING 

We rewrite the expression for the 4>j in terms of the non-spherical part of local density: 
ij 1 - 1 ™ 3 = 6 q {^iY cos 9 + 4>-j\ x sin 6 cos ip + 4>-p t v sin 6 sin ip^j 
~ —B^^R J' dr'R, r Ke ns (r' \t) - ^-D(T,0,<f>) 

-Ib^(0°)-L jf di J R, r &KE na {r',t)-lF(4 )-^ J dr'R, r R 3 Ke ns (r',t) 

cos 8 + D x sin 9 cos cp + D y sin^sintp. By examining the dilaton equations of motion in the 
there is a component of the leading-order (8, (^)-dependent term in the exterior expansion or 

1 [ r 

--B^((j) c )R / dr'R, r K,e ns (r',t) 

Joe 



where D(T,6,(f>) := D z cos 
FRW region, we can see 
<j) behaves like 
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for R <C H 1 and t fixed. Therefore matching requires that we choose D such that 

5i((> { p ns = 5 q ^ t z cos9 + (p { p x sm9cosip + <f>P t v sm9sm^ (17) 
~ ~B >4> (</><}) R [ dr'R r Ke ns {r',t) + \b, 0(^)^2 / dr' R, r R 3 K e ns (r' , t) (18) 

Joo " JR=R S 

-\ F {^)^j dr'R. r R 3 K e ns (r',t). 

n JR=R* 



IR=R S 

The interior expansion is now fully specified to order 0(5 P ). We are interested in the behaviour of <pi it and we find 



(i)n« 2 
j,t ~ 3 



i?,, (0?) i? ^ dr'R, r R, t K£ns ^ t] + \b i4 , (0?) f dr , i2, I .iJ,ti2 2 Ke na (r' ) t) (19) 

Joo it 6 K J R=Rs 

~F (4>o) ^ / dr'R, r R, t R 3 K e ns (r', t) + (4> ) RR )t Ke ns (r, t). 

This expression is valid whenever R s ^> 2m, and the requirements for matching are satisfied. In these cases we expect 
F (cf>o) w B t ((/>5) + 0(2m/R s ); so, approximately, we have 

<^" S ~ ^ (0c) i? f dr'^fl ^"'^^ + \b i4 , (4> c ) RR, t Ke ns (r, t). 

6 Joe H 6 

In the case, where R s = 2m, and our 'star' is actually a black-hole, we require F (0o) to ensure that the <f> is well- 
defined as R — ► 2m. Even so, in this case, equation (16) will not be strictly valid, since it was derived under the 
assumption of R s » 2m. By inspection of the dilaton evolution equation in the interior, eq. (14), however, we expect 

that 5 q (j)"p t ns near the black-hole horizon to be of similar magnitude to the RHS of eq. (16). 
Combining the results of this paper with those for the spherically symmetric case we find: 



'>i,t ~ cp c ,t 



B^(<pc) f dr'R ir R it K^e s {r',t) + \B^{(l) c )R [ dr'R r R t ' ^ +\B t4> (<j> c )RR tt Ke na (r,t). 



We require that \{4>i,t — <Pct) / <Pc,t \ <C 1 for 1 to hold and so ensure that local observations will detect variations of < 
occurring on cosmological scales. 



VII. GENERALISATION: A CONJECTURE 



So far, we have found an analytic approximation to the values of <fi and <fi c j in the interior. More succinctly (although 
less explicitly) we can say that, to leading order in 5, the values of <fi, <j> t t and (f>^ r can all be found everywhere outside 
the exterior region from the approximation: 

<j)~<p h om{t) + <Pl{$,t), (20) 

where (pi is the solution to: 

O sch (pi = B^nAe 

with Ae = e(x,t) — e c (t), O sc h is the wave operator in a Schwarzschild background, and t is the proper time of a 
comoving observer. This is solved w.r.t. the boundary conditions (pi — > as R — > 00 (where R = is the centre of 
our 'star') and the flux out of the 'star' is as given by equations (12) and (13). The homogeneous term is 



<t>hom{t) =4> c (t + At(x, a;)) 

where the lag, At(x,t), is defined by: 

V 2 At - v* ■ V(v* • VAi) - v* • VAt(V • v*) = -V • Av, 

with V, — di, i — {1, 2, 3}, and Aw = v — Hx, where v is the velocity of the dust particles relative to R = \\x\\ = 0. 
The velocity v* has the following properties: v* — v in some region that includes all the interior and excludes all of 
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the exterior; v* — Av everywhere else. In a general sense, the interior and exterior are two disjoint regions of total 
spacetimes where general-relativistic effects are non-negligible at leading order (e.g. such as when w 1). The 
interior region should be closed, and in the exterior region ||Au|| is small. So, v* should be defined in such a way that 
it respects all the symmetries of the spacetime and so that 1 1 -tT 5 " s * I ! "C 1 everywhere outside the interior region. This is 
required to ensure that At, as defined above, is finite. It can be seen to come out of the matching procedure. When 
the background spacetime satisfies the conditions given below, the precise way in which v* is defined does not effect 
the leading order behavior of At. For boundary conditions, we must require the flux out of At out of the 'star' to 
vanish, and require At — > as Av — ► 0, i.e. as R — ► oo. This is the natural generalisation of what has been seen in the 
Szekeres-Szafron backgrounds v 1 — R.t(R,t)5 l R . In these cases the equation is just an ordinary differential equation 
in R with solution: 

f A p/ (R,t (Rf, t) - HR' + R. t (R s , t) + HR a ) f°° , (R, t (R 1 , t) — HR' + R, t (R s , t) + HR a ) 



(R',t) J A 1 - (R it (R',t) -HR) 2 

where A is some arbitrary value of A in the intermediate region, and each A represents a particular choice of definition 
for v*. This expression is only valid to leading order in the interior and intermediate regions. To this order all 
choices for A are equivalent. Near R = R s , to leading order in S = Lj/Le, this ensures that d(f + At) ~ dv, 
where v — t sc h + R + 2m \n{R/2m — 1) is the advanced time coordinate and t sc h is the standard, curvature-defined, 
Schwarzschild time-coordinate. The solution for (frhom is then, to leading order in 5, just the particular one given by 
Jacobson in [29]. We have assumed that the generalisations of the Szekeres-Szafron result for 4> hold. We have only 
proved that this assumption holds for the subset of Szekeres-Szafron spacetimes for which the matching procedure 
works. Nonetheless, based on this analysis, we conjecture that 20 provides a good numerical approximation to the 
value of (p, and by differentiating once, to 0,t and di<j>, i — {1, 2, 3}, near the surface of our 'star', for any dust plus A 
spacetime that can be everywhere considered to be a weak-field perturbation of either Schwarzschild, Minkowski, or 
FRW spacetime; that is, 

R 2 nAe(R,t) < 1, 2(M(R, t) - m)/R < 1 

where M(R,t) is the gravitational mass contained inside the surface (R,t) = const. One could seek to motivate our 
conjecture as some sort of analytical continuation from the Szekeres-Szafron spacetimes to more general backgrounds, 
but such arguments would, we believe, be hard to frame in any rigorous context and are beyond the scope of the 
analysis in this paper. 



VIII. DISCUSSION 



In this paper we have extended the analysis of [13] to include a class of dust-filled spacetimes without any symmetries 
provided by the Szekeres-Szafron metrics. Again, we have used the method of matched asymptotic expansions to link 
the evolution of the dilaton field, <p, in an approximately Schwarzschild region of spacetime to its evolution in the 
cosmological background. By these methods, we have provided a rigorous construction of what has been simply 
assumed about the matching procedures in earlier studies [30]. We have also analysed, more fully, the conditions that 
we need the background spacetime to satisfy for the matching procedure to be valid, and we have interpreted these 
conditions in terms of their requirements on the local energy density. Finally, we have conjectured a generalisation of 
our result to more general spacetime backgrounds than those considered here. 

By combining the results found here with those of the previous paper, we conclude that, in the class of quasi- 
spherical Szekeres spacetimes in which the matching procedure is valid, the local time variation of the dilaton field 
will track its cosmological value whenever: 

(ge) C dr'R, rK A(R, t e a (r', t)) + f £ (&) R £ Ar' R. r R. t gjgj + \B± (</> c ) RR, t ne ns (r, t) 

<i>c 

When the cosmological evolution of (f> is dominated by its matter coupling: 4> c <~ 0(B ^H^ 1 ks c ), this condition is 
equivalent to: 

Ho r^ r A(W..» + 1 r, - l£=i (£i) + U, R /^A « i. 

Joo £cit) 3 s c (t) 3 e c (t) 

In the other extreme, when the potential term dominates the cosmic dilaton evolution, the left-hand side of the 
above condition is further suppressed by a factor of B^(cj) c )/V^((j) c ) -C 1. As in our previous paper, [13], we can see 



< 1. (21) 
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that for a given evolution of the background matter density, condition (1) is more likely to hold (or will hold more 
strongly) when \B^((j) c )KS c /V t( p(4> c )\ <C 1. We reiterate our previous statement that: domination by the potential term 
in the cosmic evolution of the dilaton has a homogenising effect on the time variation of tfi. 

The non-spherically symmetric parts of energy density enter into the expression differently. The magnitude of 
the terms on the left-hand side of eq. (21) is, as in the spherically symmetric case, still (H RAR^e/e) (R,t) where 
(•) (R,t) represents some 'average' over the region outside the surface (R,t) = const. We should note that, given the 
condition on ne that has been required for matching, the leading-order contribution to ne ns is everywhere of dipolc 
form and this is responsible for the special form of the average over the non-spherically symmetric terms. We can also 
see that, as a result of form of eqn. (21), peaks in ne ns that occur outside of the interior region will, in the interior, 
produce a weaker contribution to the left-hand side of eqn. (21) than a peak of similar amplitude in a spherically 
symmetric energy density ne s . This behaviour would continue if we were also to account for higher multipole terms 
in K£ ns . The higher the multipole, the more 'massive' the mode, and the faster it dissipates. 

If we are interested in finding a sufficient condition (as opposed to a necessary and sufficient one) for (1) to hold 
locally, then in most circumstances we will be justified in averaging over the non-spherically symmetric modes in 
the same way as we average over the spherically symmetric ones. In most cases, this will over-estimate rather than 
under-estimate the magnitude of the left-hand side of our condition, (21). This reasoning leads us to the statement 
that for <p(x,t) w <j> c (t) to hold locally it is sufficient that: 



where dl := dri? r , v = R :t is the velocity of the dust particles, lim^^oo v = H R. We make the same generalisation 
that we did in Paper I by taking j(R) to run from R to spatial infinity along a past, radially-directed light-ray. In this 
way, we incorporate the limitations imposed by causality. We should also assume that the above expression includes 
some sort of average over angular directions; to be safe we could replace e by its maximum value for fixed R and t. 
This sufficient condition, (22), is precisely the generalised condition proposed in our first paper on this issue. The 
inclusion of deviations from spherical symmetry, therefore, has little effect of the qualitative nature of the conclusions 
that were found in [13]. If anything, we have seen that the non-spherical modes dissipate faster and, as a result, will 
produce smaller than otherwise expected deviations in the local time derivative of <f> from the cosmological ones. 

On Earth we should expect, as before, that the leading-order deviation of <j> from <j> c is produced by the galaxy 
cluster in which we sit, and that for a dilaton evolution that is dominated by its coupling to matter, this effect gives 
I~6x lO -3 ^" 1 /! <C 1, where Q m w 0.27 and h w 0.71. If the cosmic dilaton evolution is potential dominated then 
1 is even smaller. We conclude, as before, that irrespective of the value of the dilaton-to-matter coupling, and what 
dominates the cosmic dilaton evolution, that 



will hold in the solar system in general, and on Earth in particular, to a precision determined by our calculable 
constant X. We also conclude, as before, that whenever 2" <C 1 near the horizon of a black hole, there will be no 
significant gravitational memory effect for physically reasonable values of the parameters [29, 31]. 

Our result relies on one major assumption: the physically realistic condition that the scalar field should be weakly 
coupled to matter and gravity - in effect, the variations of 'constants' on large scales must occur more slowly than 
the universe is expanding and so their dynamics have a negligible back-reaction on the cosmological background 
metric. In this paper we have removed the previous condition of spherical symmetry at least in as far as the spacetime 
background is well described by Szekeres-Szafron solution. We have therefore extended the domain of applicability 
our general proof: that terrestrial and solar system based observations can legitimately be used to constrain the 
cosmological time variation of supposed 'constants' of Nature and other light scalar fields. 




(22) 



0(x,t) « Mt) 
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